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Abstract 


Box-counting  and  EDM  methods  were  used  to  measure  the  fractal  dimension  of  round  turbulent  cryogenic  nitrogen  jets  at  pres¬ 
sures  ranging  from  subcritical  to  supercritical  pressures.  Both  method  produced  similar  trends,  with  close  quantitative  agree¬ 
ment  for  a  suitably  small  box-counting  scale.  At  subcritical  pressures,  the  fractal  dimension  was  found  to  be  consistent  with  the 
fractal  dimension  of  a  spray  in  the  2nd  wind-induced  atomization  regime.  The  fractal  dimension  tended  to  increase  as  pressure 
increased,  until  at  supercritical  pressures  the  fractal  dimension  was  found  to  be  consistent  with  that  of  gas  jets  and  mixing  lay¬ 
ers.  The  results  constitute  additional  quantitative  evidence  for  the  hypothesis  that  subcritical  jets  exhibit  mainly  spray-like  be¬ 
havior,  while  supercritical  jets  exhibit  mainly  gas-like  behavior.  This  appears  to  have  been  the  first  time  pressure  effects  on  the 
fractal  dimension  of  turbulent  jets  has  been  measured. 


Introduction 


As  combustion  chamber  pressures  become  increasingly  higher  in  order  to  realize  performance  and/or  efficiency  benefits  in  a 
wide  range  of  propulsion  and  energy  conversion  applications,  the  injected  fluid(s)  may  experience  ambient  pressures  exceeding 
the  critical  pressure(s)  of  the  propellants.  For  example,  in  the  cryogenic  liquid  hydrogen/liquid  oxygen  Space  Shuttle  Main 
Engine,  the  thrust  chamber  pressure  is  about  22.3  MPa,  and  thrust  chamber  pressures  for  the  oxygen/hydrogen  Vulcain  engine 
used  to  launch  the  Ariane  5  vehicle  have  been  recorded  to  reach  up  to  10  MPa.  Both  of  these  pressures  exceed  the  critical  pres¬ 
sure  of  5.043  MPa  for  liquid  oxygen.  In  many  applications,  including  liquid  rockets,  the  initial  temperature  of  the  propellant  is 
often  below  the  critical  temperature,  and  later  heats  up  to  a  supercritical  temperature  as  it  is  mixed  and  burned  in  the  combus¬ 
tion  chamber. 

As  chamber  pressures  approach  and  exceed  the  critical  pressure,  combustion  mechanisms  can  be  expected  to  become  distinctly 
different  than  conventional  spray  combustion  processes  at  low  pressures.  For  example,  for  single  component  fluids,  the  distinct 
difference  between  the  gas  and  liquid  phases  disappears  when  the  pressure  exceeds  the  critical  pressure  or  the  temperature  ex¬ 
ceeds  the  critical  temperature.  Surface  tension  and  the  enthalpy  of  vaporization  vanish,  and  large  variations  in  the  density, 
thermal  conductivity,  and  mass  diffusivity  occur  near  the  critical  point.  For  multicomponent  fluids,  the  solubility  of  the  gas 
phase  in  the  liquid  phase  increases  as  pressure  approaches  the  critical  pressure,  and  mixture  effects  need  to  be  taken  into  ac¬ 
count  in  calculating  the  critical  properties.  When  mixture  effects  are  considered,  the  critical  pressure  of  a  mixture  can  be  found 
to  be  several  times  the  critical  pressure  of  the  propellant  of  interest.  Understanding  and  modeling  the  effect  of  these  differences 
on  the  shear  layer  of  the  injected  propellant  are  vital  to  understanding  the  mixing  and  performance  characteristics  of  these  high 
performance  engines. 

Until  recently,  understanding  of  the  injection  process  under  these  conditions  has  been  limited.  However,  recent  experimental 
and  theoretical  work,  motivated  to  a  large  extent  by  liquid  rocket  combustion  applications,  has  begun  to  allow  a  picture  of  the 
shear  layer  dynamics  to  emerge.  See  for  example,  Mayer  et  al  [1,2],  Chen  and  Sui  [3],  Woodward  and  Talley  [4],  Harstad 
and  Bellan  [5],  Delplanque  and  Sirignano  [6],  Oefelein  and  Yang  [7],  and  Chehroudi  et  al.  [8,  9].  One  of  the  important  ex¬ 
perimental  findings  was  the  transition  in  the  morphology  of  the  injected  jet  boundary  as  chamber  pressure  increases  from  a 
subcritical  to  supercritical  pressure.  Chehroudi  et  al  [8,9]  studied  cryogenic  liquid  nitrogen  jets  (simulating  liquid  oxygen  in 
liquid  rocket  engines)  injected  into  room  temperature  nitrogen  at  various  subcritical  and  supercritical  pressures.  The  initial 
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.  temperature  of  the  jet  was  between  90  and  110  K,  which  is  below  the  critical  temperature  of  126.2  K  for  nitrogen.  Back- 
illuminated  visualizations  of  the  results  are  shown  in  Fig.  1.  It  was  found  that  at  subcritical  chamber  pressures,  the  jet  appears 
to  be  liquid-like  with  instabilities  that  grow  downstream  of  the  injector.  In  the  shear  region,  very  fine  drops  are  found  to  be 
ejected  from  the  jet.  Major  structural  changes  occur  as  the  pressure  is  increased  to  near  the  critical  pressure.  Drops  are  no 
longer  detected  in  the  shear  region.  These  drops  are  replaced  by  small  comb-like  ligaments  which  “dissolve”  into  the  surround¬ 
ing  media  due  to  the  combined  effects  of  the  reduction  in  surface  tension  and  negligible  latent  heat  of  vaporization.  As  cham¬ 
ber  pressure  is  increased  further  beyond  the  critical  pressure,  the  jet  begins  to  take  on  the  appearance  of  a  turbulent  gas  jet  in¬ 
jected  into  a  gaseous  environment.  A  quantitative  examination  of  the  shear  layer  growth  rate,  using  both  high-speed  back-lit 
images  and  Raman  scattering  data,  also  indicates  that  the  growth  rate  of  the  jet  was  consistent  with  that  of  round  sprays  at  sub- 
critical  pressures,  and  consistent  with  that  of  shear  layers  surrounding  gas  jets  at  supercritical  pressures.  See  Chehroudi  et  al 
[10]  for  more  details. 

This  metamorphic  change  of  the  cryogenic  jet  or  shear  layer  boundary  was  the  motivation  for  investigating  the  change  in  the 
fractal  dimension  of  the  boundary  in  the  present  study.  The  objective  was  to  further  quantify  the  extent  to  which  the  jets  may 
behave  like  conventional  sprays  at  subcritical  pressures,  but  like  gas  jets  at  supercritical  pressures.  The  present  study  appears  to 
be  the  first  time  the  fractal  dimension  of  jets  has  been  investigated  over  a  range  of  subcritical  to  supercritical  pressures. 


The  Fractal  Dimension 

The  notion  of  a  fractal  has  been  given  a  strong  foundation  and  application  by  its  founder  Mandelbrot  [11],  and  is  intimately 
connected  to  the  concept  of  scaling  and  self-similarity.  Peitgen ,  et.  al  [12]  pointed  out  that  in  examining  the  structure  of  a  cau¬ 
liflower,  the  head  may  be  found  to  be  composed  of  several  branches  or  parts,  which  when  removed  and  compared  with  the 
whole  are  very  much  similar  in  shape  to  the  whole,  only  smaller.  The  parts  can  in  turn  be  decomposed  into  even  smaller  parts, 
which  again  are  very  much  similar  in  shape  to  the  whole  as  well  as  to  the  first  generation  of  parts.  This  self-similarity  as  a 
function  of  scale  only  carries  through  for  about  three  to  four  stages,  after  which  the  cauliflower  is  no  longer  self  similar.  The 
scales  over  which  self  similarity  applies  is  said  to  range  from  an  inner,  or  smallest,  scale,  to  an  outer,  or  largest,  scale. 

The  concept  of  dimension  in  mathematics  has  a  historical  evolution.  In  Euclidean  geometry,  a  line,  a  plane,  and  a  surface  in 
space  have  dimensions  of  one,  two,  and  three  respectively.  In  topology,  things  are  slightly  different.  For  example,  a  straight 
line  and  any  non-self-intersecting  curve  are  equivalent  and  have  the  same  “topological  dimension.”  Finally,  these  concepts 
were  shaken  by  the  introduction  of  the  so-called  “space-filling”  curves.  Given  some  patch  of  a  plane,  there  exists  a  “space¬ 
filling”  curve  which  meets  every  point  in  that  patch.  Thus,  a  curve,  which  by  nature  is  something  one-dimensional,  can  fill 
something  two-dimensional.  This  questions  the  intuitive  perception  of  curves  as  one-dimensional  objects,  because  they  fill  the 
plane,  an  object  which  is  intuitively  perceived  as  two-dimensional.  Looking  at  a  self-similar  curve  with  the  notion  defined  in 
relation  to  the  cauliflower,  there  is  a  power  law  relation  between  the  number  of  pieces,  N,  and  the  reduction  factor,  e,  as: 
N  =  (l/£)°s  .  In  other  words,  this  power-law  relationship  states  that  if  a  self- similar  curve  is  reduced  by,  say,  8  =1/2,  the 

number  of  similar  pieces  increases  by  2Ds  .  The  quantity  D§  is  defined  as  the  “self- similarity  dimension”  of  the  curve.  Obvi¬ 
ously,  for  a  line,  a  square,  and  a  cube,  we  have  Ds=l,  2,  3,  respectively,  which  happen  to  be  the  topological  dimensions.  For 
self-similar  curves  this  dimension  is  determined  through  an  exercise  to  measure  the  total  length  of  the  curve  using  rulers  with 
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different  compass  settings  (or  yardstick  lengths).  The  reciprocal  compass  setting  can  be  considered  to  be  the  precision  of  the 
length  measurement.  Mechanistically,  one  selects  a  fixed  compass  setting  and  attempts  to  find  the  number  of  compass  settings 
needed  to  cover  the  entire  curve.  In  other  words,  one  measures  the  total  length  of  the  curve  using  a  fixed  setting  of  the  com¬ 
pass.  This  is  repeated  for  different  compass  settings.  Plotting  the  number  of  compass  lengths  required  to  cover  the  curve  at  a 
given  compass  setting  as  a  function  of  the  reciprocal  compass  setting  itself  traces  a  line  in  a  log-log  plot  for  a  self-similar 
curve.  The  slope  of  this  line,  d,  and  the  self-similarity  dimension,  Ds,  can  be  shown  related  by:  Ds=  1  +  d,  where  d  is  positive 
in  magnitude.  Usually,  however,  the  number  of  compass  lengths  at  a  given  fixed  setting  is  plotted  against  the  compass  setting 
itself  (which  in  this  case,  Ds=  1  -  d  as  d  is  negative  in  magnitude).  Motivated  by  this  result  for  self-similar  curves,  the  fractal 
(compass)  dimension  for  shapes  that  are  not  strictly  self-similar  curves  such  as  the  jet  boundary  we  would  like  to  analyze  are 
defined  as:  Dc=  1  +  d.  With  this,  the  fractal  dimension  of  the  coast  of  the  England  that  has  many  wrinkles  is  1.36  whereas  the 
that  of  the  state  border  of  Utah  having  smooth  straight  lines  is  equal  to  1 .  Note  that  a  non-integer  number  is  attributed  to  the 
dimension  (fractal).  The  fractal  dimension  of  any  curve  is  between  1  and  2  and  the  more  wrinkled  and  “space-filling”  it  looks 
the  larger  the  dimension  value  becomes.  Natural  curves,  similar  to  a  cauliflower,  are  self-similar  only  to  within  a  narrow  range 
of  scales.  Our  objective  here  is  to  measure  the  fractal  dimension  of  the  interface  of  the  jets  injected  into  the  chamber  to  exam¬ 
ine  if  any  pattern  is  uncovered.  To  our  knowledge,  this  is  the  first  application  of  the  fractal  approach  to  the  liquid  jet  interface 
under  high  pressures  up  to  supercritical  conditions. 

Practical  Measurement  of  the  Fractal  Dimension 

In  the  literature  there  are  several  methods  to  determine  the  fractal  dimension  of  a  boundary.  For  example  the  above-mentioned 
ruler  method  is  also  called  the  caliper  method  or  the  arc-swinging  method.  This  method  is  difficult  to  implement  for  computer 
analysis,  so  more  computer-compatible  and  faster  algorithms  have  been  devised.  One  popular  method  for  evaluating  the  fractal 
dimension  is  referred  to  as  box-counting,  or  the  Kolmogrov  dimension  (see  Russ  [13],  [14],  and  Kaye  [15]).  One  overlays  the 
curve  on  a  square  grid  with  mesh  size  £  and  counts  the  number  of  grid  boxes  which  contain  any  part  of  the  curve.  This  gives  a 
number  N(e)  which  depends  on  the  mesh  size  e.  The  negative  of  the  slope  of  the  line  fitted  to  the  linear  section  of  the  plot  of 
log[N(e)]  vi-.  log(£)  is  the  box-counting  dimension.  Another  approach  leads  to  what  is  referred  to  as  Minkowski  dimension,  see 
Russ  [13].  In  this  method,  the  center  of  a  circle  with  radius  £  is  swept  continuously  along  the  curve.  This  creates  what  is 
known  as  a  Minkowski  sausage  around  the  curve.  The  area  of  the  sausage  as  a  function  of  the  radius  of  the  circle  is  drawn  on  a 
log-log  plot.  The  Minkowski  dimension  is  equal  to  (2  -  the  slope  of  the  linear  region).  In  this  paper  both  box-counting  and 
Minkowski  methods  are  used  to  calculate  the  fractal  dimension.  For  the  latter  method,  the  Euclidean  Distance  Mapping  (EDM) 
algorithm  of  Russ  [13]  is  implemented.  Consider  a  black  object  in  a  white  background  where  the  fractal  dimension  of  the 
boundary  between  the  two  is  desired.  EDM  is  an  image  processing  operation  that  is  applied  to  a  black  and  white  image  to  pro¬ 
duce  a  gray-scale  image  in  which  each  black  pixel  is  given  a  brightness  value  equal  to  its  distance  from  the  nearest  white  back¬ 
ground  point.  Therefore,  the  larger  the  distance  of  a  pixel  from  the  nearest  border  the  higher  its  assigned  gray-scale  intensity 
value,  being  linearly  proportional  to  the  distance.  The  same  procedure  can  assign  values  to  the  white  background  pixels  based 
on  their  distance  from  the  nearest  point  on  the  object.  Thresholding  the  distance  map  for  either  the  object  or  the  background 
produces  uniform  erosion  and  dilation  of  the  boundary  to  any  desired  distance  from  the  original  boundary,  without  iteration. 
The  EDM  method  is  superior  to  the  classical  erosion/dilation  method,  see  Russ  [14].  It  is  really  not  necessary  to  perform  the 
thresholding  operations.  Counting  the  pixels  as  a  function  of  their  brightness  (distance)  values  produces  a  plot  that  directly 
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provides  the  Minkowski  dimension.  Both  Box-counting  and  Minkowski  methods  are  shown  to  provide  reliable  results,  see  Hall 
et  al  [16].  However,  the  EDM  algorithm  in  determining  the  fractal  dimension  is  more  isotropic,  more  efficient,  and  generally 
more  precise  than  the  box-counting  method,  and  easy  to  implement  on  a  computer,  see  Russ  and  Russ  [17].  Both  techniques 

are  limited  to  only  self-similar  and  not  self-affine  boundaries. 

/ 

Fractal  Analysis  of  Cryogenic  Jets 

Before  running  the  fractal  analysis  program  for  the  computation  of  fractal  dimensions,  the  images  need  to  be  conditioned.  In 
particular,  corrections  to  images  are  applied  to  account  for  the  reference  background  condition.  Then  smoothing  and  twice- 
sharpening  filters  are  applied.  Also,  an  appropriate  pixel  intensity  threshold  level  (described  later)  is  selected  to  identify  the 
boundary  of  the  jet  injected  into  the  chamber.  A  200X200  pixel  image  of  a  square  boundary  within  our  432  x  477  full  frame 
size  is  used  as  a  test  image  to  evaluate  the  fractal  analysis  program.  Theoretically,  the  fractal  dimension  of  the  smooth  bound¬ 
ary  of  this  square  must  be  equal  to  one.  The  analysis,  however,  gives  a  fractal  dimension  higher  than  one  by  at  most  3%  for  the 
box-counting  and  EDM  methods.  For  the  box-counting  dimension,  two  different  methods  are  used  to  change  the  mesh  size.  In 
one  case,  the  mesh  size  is  increased  by  the  following  sequence,  lxl,  2x2,  3x3,...,  32x32  pixels,  and  in  the  other  by  2x2,  4x4, 
6x6,  8x8,  64x64  pixels.  In  this  paper  they  are  referred  to  as  BOX32  and  BOX64  respectively. 

In  order  to  define  a  closed-loop  jet  boundary  for  fractal  analysis  of  the  jet  images  in  Fig.  1,  one  needs  to  manually  select  a  sec¬ 
tion  of  the  jet  image  that  excludes  regions  covering  the  image  of  injector  tip  and  the  part  of  the  image  near  the  lower  edge  of 
the  frames.  This  creates  two  smooth  straight  lines  joining  the  left  and  the  right  boundaries  of  the  jet  at  the  top  and  the  bottom. 
The  effect  of  adding  these  two  smooth  sections  on  the  fractal  dimension  value  is  analyzed  by  several  different  methods.  First, 
the  left  and  the  right  boundaries  are  manually  moved  close  to  each  other  so  the  length  of  the  upper  smooth  line  is  minimized. 
The  minimization  is  done  on  the  upper  section  because  the  jet  diverges  downward  in  the  images.  Second,  the  upper  and  lower 
smooth  lines  are  replaced  by  pieces  of  the  right  and/or  left  boundaries  of  the  jet.  The  fractal  dimension  increased  by  2  to  3  % 
applying  such  corrections  on  the  jet  images  for  various  chamber  pressures.  For  ease  of  analysis,  however,  we  determine  the 
fractal  dimension  of  jet  images  with  the  two  upper  and  lower  smooth  lines  and  uniformly  apply  a  2.5%  increase  to  all  final  val¬ 
ues  as  a  first-degree  correction.  The  selection  of  the  pixel  intensity  threshold  for  jet  boundary  determination  is  an  important 
factor.  For  this  reason,  at  each  tested  chamber  pressure  the  fractal  dimensions  of  the  jets  are  calculated  for  a  wide  range  of 
threshold  levels  and  the  results  are  plotted.  A  maximum  excursion  of  up  to  10%  is  observed  at  high  pressures  as  the  threshold 
is  varied.  In  most  cases,  one  can  detect  a  plateau  region  in  which  a  threshold  value  can  be  selected.  If  not,  a  value  at  the  center 
of  the  explored  threshold  range  is  chosen.  As  a  cross  check,  the  boundary  generated  at  the  chosen  threshold  is  visually  overlaid 
and  inspected  to  follow  the  original  gray  scale  image  as  close  as  is  discernable.  Figure  2  shows  a  log-log  plot  of  the  number  of 
boxes  counted  versus  the  box  size  for  a  selected  jet  at  high  pressure  and  for  five  different  threshold  levels.  On  the  thresholded 
images  one  observes  many  vein-like  structures  that  are  not  part  of  the  closed-loop  jet  boundary  but  attached  to  it.  These  are 
removed  by  an  appropriately  selected  sequence  of  erosion/dilation  image  processing  cycles.  Special  treatment  is  also  required 
for  images  at  low  pressures  (i.e.  the  subcritical  regime).  The  boundary  of  the  jet  is  determined  such  that  all  isolated  and  single¬ 
pixel  touching  ligaments  and  drops  are  excluded.  This  is  achieved  by  application  of  series  of  erosion  and  dilation  steps  as  de¬ 
scribed  in  Russ[  13]. 
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In  the  past  ten  years  a  number  of  applications  of  fractal  analysis  have  been  demonstrated  in  different  disciplines.  For  example, 
Sreenivasan  and  Meneveau  [18]  computed  the  box-counting  and  Minkowski  fractal  dimension  of  the  turbulent  boundary  of  a 
smoke-tagged  turbulent  gas  jet  visualized  across  its  diameter  using  a  laser  sheet,  where  the  laser  sheet  thickness  was  smaller  or 
comparable  to  Kolmogorov  scale  to  minimize  “fuzzing”  due  to  the  finite  thickness  of  the  sheet.  Sreenivasan  [19]  summarized 
the  results  and  reported  fractal  dimensions  of  1.35,  1.34,  and  1.38  for  a  round  gaseous  jet,  a  plane  gaseous  mixing  layer,  and  a 
boundary  layer  flow,  respectively,  using  two-dimensional  digital  imaging  from  LIF  visualizations.  At  atmospheric  pressure, 
Shavit  and  Chigier  [20]  established  a  relationship  between  the  peak  of  the  fractal  dimension  calculated  at  different  streamwise 
directions  and  the  breakup  point  as  well  as  the  intact  length  of  a  jet  (or  spray)  produced  by  a  coaxial  air-assist  injector.  They 
report  the  peak- value  dimension  to  range  from  1.12  to  1.32  as  the  air  speed  is  increased  from  30  to  90  m/s.  It  is  also  interesting 
to  mention  that  Mantzaras  et  al  [21],  who  measured  the  Minkowski  dimension  of  flames  in  spark-ignited  engines,  and  North 
and  Santavicca  [22],  who  measured  the  Minkowski  dimension  of  atmospheric  pressure  premixed  flames,  both  report  a  pro¬ 
gressive  increase  in  the  fractal  dimension  from  1.05  to  1.36  as  the  turbulence  level  is  increased.  This  indicates  a  progressive 
wrinkling  and  space-filling  of  the  flame  front  by  turbulence.  Finally,  modeling  turbulent  flames  using  fractal  concepts  is  dis¬ 
cussed  in  Gouldin  [23]  and  Gouldin  et  al  [24]. 

Fractal  dimension  measurements  of  the  nitrogen  jet  images  are  plotted  as  a  function  of  the  reduced  pressure  (chamber  pressure 
divided  by  the  critical  pressure  of  nitrogen)  in  Fig.  3.  Separate  curves  are  shown  for  the  BOX32,  BOX64,  and  EDM  methods, 
and  a  fourth  curve  plots  the  average  of  all  three.  The  value  of  the  fractal  dimension  at  each  pressure  is  the  average  of  18  to  20 
images.  The  BOX64  method  produces  systematically  higher  measurements  than  the  other  two,  which  is  believed  to  be  due  to 
contributions  beyond  the  32x32  grid  size.  Considering  the  better  accuracy  of  the  EDM  method  and  its  closeness  to  the  BOX32 
results,  the  latter  two  methods  are  considered  to  produce  more  representative  measurements  of  the  fractal  dimension.  Also  plot¬ 
ted  in  Fig.  3  are  fractal  dimensions  measured  by  other  researchers  of  liquid  sprays  and  gaseous  jets,  mixing  layers,  and  bound¬ 
ary  layers.  The  measurements  reported  for  Dimotakis  et  al.  [25]  and  Taylor  and  Hoyt  [27]  are  our  measurements  made  from 
scanned  images  of  pictures  contained  in  those  references. 

The  general  trend  of  all  the  curves  in  Fig.  3  is  that  the  fractal  dimension  increases  as  the  chamber  pressure  increases.  At  low 
pressures,  the  trend  of  the  curves  is  towards  the  Euclidian  value  of  1  for  a  smooth  circular  cylinder  with  no  surface  irregulari¬ 
ties.  The  magnitude  of  the  measurements  at  the  lowest  pressures  where  the  jets  are  subcritical  is  consistent  with  the  fractal  di¬ 
mension  of  a  2nd  wind-induced  spray  measured  from  the  scanned  images  of  Taylor  and  Hoyt  [27].  As  pressure  increases,  a  dip 
in  the  curves  occurs  at  a  reduced  pressure  of  about  0.9,  due  to  a  transition  out  of  the  2nd  wind- induced  atomization  regime.  At 
the  highest  pressures  where  the  jets  are  supercritical,  the  magnitude  of  the  fractal  dimension  becomes  consistent  with  the  fractal 
dimensions  of  gaseous  jets  and  layers. 

These  results  support  the  hypothesis  that  subcritical  jets  exhibit  mainly  spray-like  behavior,  while  supercritical  jets  exhibit 
mainly  gas-like  behavior.  Earlier  support  for  this  hypothesis  included  qualitative  observations  of  the  jet  appearance  and  the 
quantitative  spreading  rate  measurements  of  Chehroudi,  et.  al  [8].  The  present  fractal  dimension  measurements  constitute  the 
second  quantitative  respect  and  the  third  overall  respect  in  which  the  hypothesis  is  true. 
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Summary  and  Conclusions 


The  fractal  dimensions  of  round  turbulent  cryogenic  nitrogen  jets  have  for  the  first  time  been  measured  at  various  pressures 
ranging  from  subcritical  to  supercritical.  The  jets  were  injected  at  an  initially  subcritical  temperature  into  room  temperature 
gaseous  nitrogen.  A  comparison  of  the  box-counting  and  EDM  methods  of  measuring  the  fractal  dimension  revealed  the  same 
trends  as  a  function  of  pressure  in  all  cases,  and  close  quantitative  agreement  when  the  box-counting  scale  is  suitably  small.  At 
subcritical  pressures,  the  fractal  dimension  is  consistent  with  a  conventional  spray  in  the  2nd  wind-induced  atomization  regime. 
The  fractal  dimension  tends  to  increase  as  pressure  increases,  until  at  supercritical  pressures  the  fractal  dimension  is  consistent 
with  those  of  gaseous  jets  and  layers.  These  results  constitute  further  quantitative  evidence  that  subcritical  jets  exhibit  mainly 
spray-like  behavior,  and  supercritical  jets  exhibit  mainly  gas-like  behavior. 
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Figure  1.  Back-illuminated  images  of  cryogenic  nitrogen  injected  into  ambient  nitrogen  at  a  fixed  supercritical  tempera¬ 
ture  of  300  K  but  varying  sub-  to  supercritical  pressures  (Political  =  3.39  MPa;  Tc  =  126.2  K).  Reduced  pressure,  Pr  = 
Pch/PCriucai=  0.23,  0.43, 0.83,  1.03, 1.22, 1.62,  2.03, 2.44,  2.74;  from  lower  right  to  upper  left.  Reynolds  number=  25,000 
to  75,000.  Injection  velocity:  10-15  m/s.  Fraud  number:  40,000  to  110,000.  Injectant  temperature:  99  to  120  K. 


Figure  2.  Log-log  plot  of  the  number  of  boxes  counted  as  a  function  of  the  box  size  for  four  different  thresholded  binary 
images  of  the  same  original  high-chamber-pressure  gray-scale  jet  image  indicating  effects  of  the  threshold  level  on  the 
fractal  dimension.  The  slope  of  the  least  square  fit  line  is  the  negative  of  the  box-counting  dimension.  R2  is  the  correla¬ 
tion  coefficient  of  the  linear  fits. 
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Figure  3.  Box-counting  and  Minkowski  (EDM  algorithm)  fractal  dimensions  of  the  visual  boundary  of  the  jet  as  a 
function  of  the  reduced  chamber  pressure  for  N2-into-N2  injection.  AVERAGE  (N2  into  N2)  is  the  average  of  the  three 
shown  fractal  dimension  measurements. 
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